Single Magnetic Atom on a Surface: Anisotropy Energy and Spin Density 



o 



(N 



Jing-Neng Yao^ and Chiung-Yuan Lin^ 

^Department of Electronics Engineering and Institute of Electronics, 

National Chiao Tung University, Hsinchu, Taiwan 

(Dated: February 24, 2013) 

Studying single-atom magnetic anisotropy on surfaces enables the exploration of the smallest 
magnetic storage bit that can be built. In this work, magnetic anisotropy of a single rare-earth atom 
on a surface is studied computationally for the first time. The single adatom and its substrate surface 
are chosen to be a Dysprosium (Dy) atom and a copper-nitrite surface, respectively, where single 
transition-metal magnetic atoms on the same surface were previously studied one atom at a time by 
scanning tunneling microscopes. We propose unconventional / and d subshell symmetries so that 
following the first-principles calculations, simple pictorial analyses of the spin-density distribution 
can be performed for the first time, independently for both a rare-earth atom Dy and a transition- 
metal Fe. The magnetic anisotropy energy of Dy on the surface is calculated to be a factor of five 
larger than the previous highest one, reaching a record-high value of 31 meV. 
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I. INTRODUCTION 

The desire of higher-density magnetic data storage is 
becoming more emergent to conform to the explosive 
growth of today's information industry. The continued 
downscahng of storage bits wih require significant en- 
hancement of magnetic anisotropy energy (MAE) per 
atom to overcome the fundamental, superparamagnetic 
limit. Single- atom magnetic anisotropy on surfaces serves 
as an important bottom-up approach of such enhance- 
ments. A surface magnetic adatom no longer has the 
isotropic orbital electronic structures of its free atom, and 
consequently its spin-orbit couplings (SOC) vary in dif- 
ferent spin orientations. Such an anisotropy is not only 
technologically relevant to the smallest magnetic stor- 
age bit that can be built, but also of great scientific 
interest for its fascinating quantum effect s^^. Pioneer- 
ing experiments^"— have been able to measure the MAE 
of single transition-metal atoms and dimers on surfaces, 
followed by computational studies^^^^^^^ that are consis- 
tent with the measurements. Among those experiments, 
scanning tunneling microscopes (STM) can even reach 
the precision such that magnetic atoms can be built, ma- 
nipulated, and measured one atom at a time at desired 
atomic sites of a CuN surface^—. There is even a very 
recent STM study^^ that can fabricate a bistable atomic- 
scale antiferromagnet, which enables a low-temperature 
demonstration of dense nonvolatile storage of informa- 
tion. However, the magnetic properties of such adatoms, 
e.g. spin density, have not yet been understood in simple, 
atomic-scale microscopic pictures. Also, MAE of non- 
transition- metal atoms, the rare earth for example, have 
not been explored, and are potential candidates with even 
higher single-atom MAE so that one may build atomic- 
scale storage bits in the future. 

We start this study from first-principles calculations 
of a single rare-earth atom Dysprosium (Dy) on a CuN 
surface. Analysis of subshell quantum numbers, orbital 
shapes, and occupations of the 4/ orbitals leads to simple 



explanation of the calculated spin-density distribution of 
the Dy atom. The success of such orbital analysis is then 
duplicated to the previously studied single Fe atom on 
the same surface^i^iii, showing that this simple picture 
works independently for both the transition-metal d and 
rare-earth / orbitals. This work of Dy on a CuN surface 
is also the first study of the MAE of a single rare-earth 
atom on a surface, which is calculated by first principles 
to be a factor of five larger than the previously largest 
single-atom MAE, Co on the Pt surface^. 



II. COMPUTATIONAL APPROACHES 

We construct a supercell of 5-layer Cu(lOO) slabs, place 
N and Dy atoms on the surface in the same way as 
the density-functional-theory (DFT) studies of single Fe 
and Mn adatoms^^^^^^, and then perform DFT calcula- 
tions using the basis set of all-electron full-potential lin- 
earized augmented plane wave^^. A naive local-density 
or generalized-gradient approximation makes the 4/ or- 
bitals be pinned at the fermi level and causes a fractional 
4/ occupation, inconsistent with the photoemission mea- 
surements of rare-earth nitride bulks, and this problem 
was shown to be resolved by performing a DFT+U cal- 
culations instead^^. To find the correct U and J values 
for DFT+U, we perform PBE+U calculation a^^i-*^^" — and 
inspect the 4/ occupation number n^f of Dy on the CuN 
surface with a series of U and J. When J = 0, and U 
ranges from to 9 eV (noting that U = denotes PBE 
itself), the occupation number stays at n^f = 9.6 regard- 
less of the U value used. This means that with solely 
U turned on, the unphysical fermi-level pinning of one 
of the / states remains the same as the naive PBE. If 
J is turned on to be a positive value 1.2 eV at /7 = 9 
eV, the / state originally pinned at fermi-level lowers in 
energy, and becomes fully occupied. However, the oc- 
cupation number then becomes n^f ~ 10, inconsistent 
with the typical n^f ?^ 9 in either a DFT+U calcula- 
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Figure 1: (a) Electron density contour of a single Dy on the 
CuN surface along the N-Dy-N raw and the out-of-plane di- 
rection, in comparison with the Fe caso^. The numbers inside 
the circles indicate the net charge on selected atoms, (b) Level 
diagrams in scale showing the calculated MAE of Dy and Fe 
on the CuN surface, respectively. The Fe case is done by 
Shick et al.^ (c) Calculated spin-density isosurfaces (second 
row, blue) of Dy on CuN at the magnitude of 0.05e/ao within 
a 3 X 3 X 3A^ cube centered at the Dy nucleus, by looking (left 
to right) along x, y, and z directions. Each stick-ball struc- 
ture corresponds to one of the three directions. The purple 
surfaces are the corresponding Fe spin density for comparison. 



tion of a DyN bulki^ or a recent quantum-chemistry, 
complete active space self-consistent field (CASSCF) cal- 
culation of a DyCl3 molecule^^. If we first determine 
[/ = 6 eV from a constraint PBE calculation^^' ^^, and 
are gradually turning on a negative J, we initially find 
that the DFT self-consistent cycles do not converge be- 
fore J reaches —0.8 eV, indicating unstable electronic 
structures. When J reaches —0.8 eV, the self-consistent 
cycles converge nicely, and give a reasonable occupation 
number 77.4/ =8.9, consistent with the typical 4/^ config- 
uration of Dy. Therefore, for Dy on the CuN surface, we 



have determined that [/ = 6 eV and J = — 0.8 eV yields 
a 4/^ configuration that agrees with both the DFT+U 
calculation of a DyN bulk^^ and the CASSCF calcula- 
tion of a DyCl3 molecule^^. Atomic charges and spins 
are calculated by Bader analysisi^S. 



III. RESULTS 

As already pointed out in the previous Fe study^ , when 
an adatom is deposited onto the Cu site of the sur- 
face, it establishes polar covalent bonds with the nearest- 
neighbored N atoms that replaces the original CuN bind- 
ing network. The calculated electron density of a Dy 
atom in the CuN surface is shown in Fig. [1^, together 
with the previously calculated Fe re-presented. As one 
can see, the Dy atom, sitting even higher on top of the 
surface, attracts its neighboring N atoms further out of 
the surface than the Fe case. We have also calculated that 
Dy and its neighboring N are +1.3 and —1.2 charged re- 
spectively. Compared with the +0.9 and —1.4 charged 
Fe and its neighboring N respectively, the Dy-N bond of 
the Dy system has a polarity approximately the same as 
the Fe-N. 

By pointing the Dy spin in the hollow, N-row, and 
out-of-plane three symmetry directions (to be called x, 
?/, and z respectively) in our PBE+U total-energy calcu- 
lations with SOC included, we obtain the Dy MAE E{n) 
of n = e^;, e^y, and e^, and compare it with the Fe case, as 
shown in the level diagram in Fig.[T|3. One notices that in 
contrast to the Fe case, the most-preferred magnetization 
axis of Dy is oriented in the out-of-plane direction, while 
both atoms have their least-preferred axis pointing in the 
hollow direction. The MAE of Dy is basically one order 
of magnitude larger than that of Fe, and is five times 
larger than the 6meV MAE of Co on the Pt surface^, the 
largest single-atom MAE reported previously. 

The previously studied Fe on the surface has 13.5% 
of spin density extends into the surrounding atoms with 
the spreading primarily along the N-row direction^. In 
contrast to Fe, when calculating Dy on the same surface, 
we find that a net spin of 6* = 2.9 is localized at the Dy 
atom, and S = 2.9 by including the spin of all atoms, 
indicating that there is no spin spreading. The S = 2.9 
of Dy includes two parts of contribution from our anal- 
ysis, S = 2.5 from the localized electronic configuration 
(4/)^ (within the Dy muffin-tin sphere), and S = 0.4 
from the delocalized {6sbd)^ (out of the muffin-tin sphere 
but within the Bader basin), where 6^5^ denotes a hy- 
bridized molecular orbital^^. The significant reduction of 
spin spreading when replacing Fe by Dy is obviously be- 
cause the Dy 4/ orbitals are more localized than the Fe 
3d. Another interesting feature is the shape of the spin 
density. When looking closely at the spin isosurfaces of 
the Dy and Fe atoms along all three crystal-symmetry di- 
rections in Fig. [It, the shapes of Dy and Fe are found to 
be approximately a hexagon and a square, respectively, 
along either the N-row or out-of-plane direction, while 
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Figure 2: Schematic plots of two different sets of 4/ orbitals, 
where blue (red) zones denote positive (negative) values, (a) 
The general self^, from left to right: (upper) yz^ ^ z{z^ — 3x^), 
x(x^ — 3y^), (lower) 2;^, z{x^—y^), y(3x^—y^). The x^z orbital 
is plotted in the inset at the right-lower corner, (b) The \Li\- 
fully-polarized set, from left to right: (upper) yz^, z{z'^ — 3x'^), 
x{x^ — 3y^), (lower) zy^, xz^, y{3x^ — y^). All six orbitals have 
a unique six-petal shape. The only exceptional shape of the 
xyz orbital is plotted in the inset at the right-lower corner. 



both atoms become more round-shape-like along the hol- 
low. This seemly mysterious observation of spin shapes 
will become clear when we later look into the Fe 3d and 
Dy 4/ orbitals in the CuN surface. 

The d orbitals in a crystal environment have the fol- 
lowing well-known subshell symmetries, 2:^, x'^ — ^^, xy^ 
yz^ zx. The exotic, complicated / orbitals, which are 
rarely shown in literatures, have two commonly used sets 
of subshell symmetries, the cubic and the general sets^^. 
As an example, we plot the 4/ general set in Fig. [2^. 
We find there are four orbitals having a unique six-petal 
shape: yz'^, z{z'^ — 3x^), x{x'^ — 3^^), and y{3x'^ — ^^), 
and are different from the rest three. The six-petal shape 
suggests that the hexagonal spin density is related to the 
4/ orbitals in some way. However, the z^ and z{x'^ — y'^) 
orbitals that do not have six-petal shapes prevent us from 



establishing a relation between the spin density and the 
general-set 4/ orbitals. Here we propose an unconven- 
tional set of / orbitals such that six of the orbitals have 
fully polarized angular momenta |I/^|, as one will see its 
advantages in the orbital analysis of the spin shapes. 
Five out of the seven orbitals, xz^, yz'^^ y{3x'^ — v'^)^ 
x{x'^ — 3?/^), and xyz still belong to the general set. The 
rest two are z{z'^ — 3x^) and zy'^, which are related to 
the general-set orbitals z^ and z{x'^ — y'^) simply by the 
following orthogonal transformation 
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where the primed coordinates are arranged as 
{x'^y'^z') = (y^z^x). It can be seen that z{3y'^ — z'^) 
and zx^, which are orthogonally transformed from the 
general-set z'^ and z{x^ — v'^)^ no longer belong to the 
general set in the unprimed coordinates, but are actu- 



ally the general-set orbitals x'{x' 
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and x'z' in 



the primed coordinates. One can easily verify that each 
orbital of this set except for xyz^ with its quantization 
axis h properly chosen along one coordinate axis x^, has 
its angular momentum fully polarized, that is, they are 
eigenstates of |L • e^l = \Li\ with eigenvalues |m| = 3 = /. 
Due to the reflection symmetry of these orbitals in the 
±6^ directions, their angular-momentum polarization has 
an sign ambiguity. The full polarization of \Li\ of these 
six 4/ orbitals can be further visualized clearly by look- 
ing at their orbital shapes. As one can see in Fig. [2]3, 
the above mentioned six fully-polarized 4/ orbitals have 
an unique six-petal shape, and can be grouped into three 
pairs. Each pair of orbitals have the same coordinate axis 
as their central symmetry axis, and are related to each 
other by exchanging their in-plane axes. The most obvi- 
ous full polarization can be realized for the y{3x'^ — y'^) 
and x{x'^ — 3y'^) orbitals, which are eigenstates of \Lz\ 
with eigenvalues \m\ = 3, and have z as their central 
axis. Similarly, each of the rest four fully-polarized or- 
bitals have their central axes as their polarization direc- 
tions. The only exception, the xyz orbital, has eight lobes 
pointing to the corners of a cube. 

Our calculations show that the Dy 4/ majority-spin 
states are all fully occupied and are very low-lying, ex- 
tremely atomic-like levels. The orbital analysis of this 
work does not depend on the details of the Dy 4/ 
majority-spin states as long as all of them are fully oc- 
cupied, and their details are not presented here. For 
minority-spin states, we plot their PDOS in Fig. [3^, and 
find that xz'^ and y{3x'^ — y'^) are mainly occupied, and 
yz'^ and x{x'^ — 3y'^) slightly occupied. The occupation 
numbers of the rest of the 4/ minority-spin states are 
negligible. The bd^s electron mentioned in a previous 
paragraph is mostly delocalized out of the Gd muffin-tin 
sphere, and has a negligible PDOS within that sphere. 
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pied. We re-calculate the PDOS and plot the minority 
spins in Fig. [Sb- The minority-spin PDOS show that 
x^ — y'^ is fully occupied, z'^ partially occupied, and the 
rest three minority-spin states basically empty. The par- 
tially occupied spin-minority z^ prohibits us from estab- 
lishing a simple picture of either the spin-density shape 
or the SOC. However, one may notice that the z^ and 
x^ — y'^ PDOS profiles depend nontrivially on the choice 
(or interchange) of the coordinate axes, while the xy^ 
yz and zx (see the insets of Fig. [3b) have trivially the 
same set of PDOS profiles with interchange of PDOS la- 
bels corresponding to the interchange of the coordinate 
axes. Therefore we search for z'^ and x'^ — y^ minority- 
spin PDOS of all three possible assignments of coordinate 
axes, i.e., the hollow, N-row, and out-of-plane directions 
are x, ?/, and z axes, respectively, and the cyclic permuta- 
tions. In all three axis assignments, all the majority-spin 
states are always fully occupied, and the xy^ yz and zx 
minority-spin always empty, while the z^ and x^ — y^ 
minority-spin states have specially simple occupations as 
shown in Fig. [3b under a particular assignment of new 
axes: the hollow, N-row, and out-of-plane directions are 
2^1, xi, and y\ axes, respectively. In this new coordinate 
system, only the z\ (or equivalently x^ in the old coor- 
dinates) orbital has paired spins, while the rest four all 
unpaired. 



IV. DISCUSSION 



Figure 3: (a) The 4/ minority- spin PDOS of Dy on the 
CuN surface without SOC. The subshells follow the choice 
of Fig. [2b. The nontrivial occupation are specified in paren- 
theses, (b) The 3(i minority-spin PDOS of Fe on the CuN 
surface. The left figure has the conventional quantization 
axes with the hollow, N-row, and out-of-plane directions be- 
ing X, ^, and z axes, respectively——. In contrast to the con- 
ventional set of axes, the right one has unconventional axes 
{xi,yi,zi) — (y,z,x), and is used to perform orbital anal- 
ysis in this work for its simple occupation numbers (either 
nearly occupied or almost empty). The PDOS of three of the 
subshells are plotted separately in each inset solely for the 
purpose that all PDOS are better visualized. 



The Dy 4/ orbitals are rather localized as indicated from 
a previous paragraph that they have an approximate 
0.4 A radius. It is a good approximation to think of 
these orbitals in the way of electronic configurations of 
an atom, and therefore the nontrivial occupation (nei- 
ther fully occupied nor completely empty) is determined 
for each orbital from the area under the curve below the 
Fermi level. The PDOS thus implies an approximately 
4/^ configuration for the Dy atom itself. 

The conventional quantization axes of the Fe 3d PDOS 
are oriented in the way that the x axis points along 
the hollow direction, y along the N row, and z out of 
plane. Such PDOS have been presented in Ref.— i^^, and 
show that all the majority-spin states are fully occu- 



With the orbital quantum numbers and occupation 
numbers determined, we now try to explain the spin- 
density shape of a Dy atom on the CuN surface. By start- 
ing with the top view, one notices that the x{x'^ — 3y'^) 
and y{3x'^ — y'^) orbitals both have the six-petal shapes 
centered about the z-axis. As we have identified from the 
PDOS analysis that y{3x'^ — y'^) has roughly paired spins, 
while x{x'^ —3y'^) has roughly an unpaired spin. The spin 
density shape along the xy plane are therefore dominated 
by the x(x^ — 3y'^) majority spin alone. Observing the 
shape of the x{x'^ — 3y'^) orbital, one then realizes that 
the hexagonal shape of the Dy spin density from a top 
view in Fig. [Th is essentially the consequence of the un- 
paired x{x'^ — 3y'^) orbital (see Fig. [4^). Similarly, the 
z{z'^ — 3x^) and xz'^ orbitals, both with a six-petal shape 
centered about the 7/-axis, contribute to the hexagonal 
spin density by looking along the y direction. However, 
the yz'^ and zy'^ orbitals, centered about the x-axis, both 
have unpaired spins. Their two six-petal shapes with 
30° relative to each other result in a slightly-round-shape 
spin-density by looking along the x direction. We can 
even further understand the unique dent at the y ends of 
the spin isosurface by the following analysis. The only 
unpaired orbital that carries a lobe in the y direction is 
yz'^^ as can be seen in Fig. [4^, and all the rest have nodes 
along y. On the other hand, both the other two axes re- 
ceive lobe contribution from two unpaired-orbit als, e.g., 
x{x^ — 3y'^) and zy'^ both contribute lobes along x. Con- 
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Figure 4: Illustration of how analyses of localized orbit als 
explain the shapes of the spin density of (a) a Dy atom on the 
CuN surface, and (b) Fe on the same surface. In both cases, 
all spin up states are occupied, and the spin-down occupation 
of each orbital is determined by the PDOS in Fig. O 



sequently, the spin isosurface shrinks in the y direction 
due to relatively less lobe contribution than the x and z 
directions. 

The simple occupation picture of the adatom's local- 
ized orbitals that explains the spin-density shape works 
not only for the rare-earth atom Dy. As we are going to 
show below, the same picture also works for the previ- 
ously studied Fe adatom. When showing the Fe PDOS 
in both the conventional quantization axes (hollow x, N- 
row y, and out-of-plane z) and the unconventional ones 
(hollow zi, N-row xi, and out-of-plane yi) in Fig. [3)3, we 
have demonstrated that in the new coordinate system 
(xi, ?/i, 2:1), only the zl orbital has paired spins, while 
the rest four are all unpaired. As seen from Fig. HJd, the 
square spin-density shapes of Fe in the CuN surface from 
the top view and N-row side view are essentially con- 



sequences of the spin-unpaired ziXi and yiZi orbitals, 
respectively, while the round shape from the hollow-site 
side view is the combination of the spin- unpaired x\ — y\ 
and xiyi orbitals. One also notice that the spin isosur- 
face has a dent only at the zi ends. Similar to the way of 
explaining the Dy spin-isosurface dent, such a dent is the 
consequence that both the xi and yi ends have lobe con- 
tributions from the x\ — y\ orbital, while zi has no lobe 
contribution mainly due to the fact that the zi orbital is 
spin-paired. 



V. CONCLUSION 

The STM moving-atom technique has demonstrated 
its capability of building, manipulating, and measur- 
ing a single atomic spin in a well-characterized environ- 
ment. First-principles calculations conclude that such an 
atomic spin forms a surface-embedded molecular mag- 
netic structure^, as well as, reproduce the measured mag- 
netic anisotropy axes^. As an ongoing study of the STM- 
engineered adatoms, this work is the first attempt to un- 
derstand the spin density of a surface magnetic atom 
in a simple, atomic-scale microscopic picture. This is 
achieved by analyzing the occupations, shapes, and an- 
gular momenta of its individual localized orbitals. These 
localized orbitals include both the well-known d orbitals 
of the transition-metal atoms and the nontrivial / of the 
rare-earth. We determine an unconventional set of 4/ 
subshell quantum numbers and an unconventional set of 
3(i quantization axes that can be used to analyze the 
shape of the spin density. The spin-density shape of a 
magnetic adatom is explained by simply counting the oc- 
cupation of each individual subshell and spin state of 
the atom's localized orbital. These studies provide an 
important microscopic picture of the atomic-scale ori- 
gins of magnetization distribution. We have also done 
the first computational study of a single rare-earth atom 
on a surface. The calculation predicts a record-high 
MAE of 31 meV. All these theoretical and computa- 
tional realizations, combined with continuing experimen- 
tal innovation^, may lead to further engineering of the 
single-atom magnetic anisotropy, further exploration of 
the smallest magnetic storage bit, and new fascinating 
applications of atomic-scale magnetism. 
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